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ABSTRACT 

A heuristic for solving large-scale zerc-one 
programming problems is provided. The heuristic is based on the 
modifications made by H. Crowder et al. (1983) to the standard 
branch-and-bound strategy. First, the initialization is modified. The 
modification is only useful if the objective function values for the 
continuous and the zero-one programming problems are close to each 
other. Given the initialization, the branch-and-bound method is 
stopped when a feasible solution to the problem is found. The 
heuristic also uses the reduced costs to fix non-basic variables to 1 
or 0. An example taken from achievement test construction illustrates 
the efficiency of the proposed heuristic. Several test construction 
problems were implemented and solved by the proposed heuristic for 
item banks with 400 items. Modifications were introduced in the LANDO 
computer program. A table illustrates that the central processing 
unit times for solving the zero-one programming problem were close to 
the times needed to solve the continuous problem. (SLD) 
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Abstract 

A heuristic for solving large-scale zero-one programming 
« 

problems is given. The heuristic is based on Crowder. Johnson 
and Padberg's (1983) modifications of the standard branch— 
and-bound strategy. First, the initialization is modified. 
The modification is only useful if the objective function 
value for the continuous and the zero-one programming problem 
are close to each other. Given the initialization, the 
branch-and-bound method is stopped when a feasible solution 
to the problem is found. The reduced costs are used to fix 
variables. As an example, the heuristic is applied to 
achievement test construction problems with good results and 
excellent CPU-times. 
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A Note on Solving Large-Scale Zero-one 
Programming Problems 

Recently, Crowder, Johnson and Padberg (1983) gave some 
methodological advances which, combined with clover branch- 
and-bound strategies, solve sparse large-scale zero-one 
programming problems in feasible time. This method may- 
replace standard branch-and-bound approaches such as in Dakin 
(1965) which for many applications take too much time 
(Papadimitriou & Steiglitz, 1982). Tn this note a heuristic 
will be given which is based on some ideas proposed by 
Crowder et al. in their branch-and-bound strategy. The 
heuristic is useful for solving large-scale zero-one linear 
programming problems with small differences between the 
coefficients in the objective function. A practical example 
from the area of achievement test construction shows that the 
heuristic solves zero— one programming problems in CPU-times 
closs to those for solving the continuous relaxations. 

Notation 

Zero-one programming problems of the following form are 
considered: 

(P) . Maximize {c'x | Ax ^ b, Xj = 0 or 1 for j = 1,.., n } 

where A is an mxn-matrix and b and c are vectors of lenghts m 
and n . 
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The continuous optimal objective function will be 
denoted by Zlp and the true lower bound on the zero-one 
optimal objective function by Z + . In the heuristic, the 
continuous optimal reduced costs . d-j . corresponding to 
variable x-j are used to fix variables at the value 0 or 1 . 

Modifications in Standard Branch-and-bound 

The modifications in the branch-and-bound strategy given 
in this section can be applied when (1) the optimal value of 
the objective function for the zero-one programming problem P 
is close to Zlp, and (2) problem P is a large-scale 
piogramming problem, it is assumed that Zlp is not equal to 
0. 

If no feasible solution to P is known, the branch-and- 
bound method is initialized by assuming Z + = -«>. But if it is 
known that the optimal value of the zero-one objective 
function is close to Zlp, the branch-and-bound method can. 
after solving the relaxation of P (< be initialized by 
Z + = K^Zlp., where is a constant (0 << <1). 

Given the above initialization it is known that every 
zero-one solution found during the search process has a value 
of the objective function between K^ZlP and Zlp- So if is 
close to 1 every solution is a good solution. This means that 
the branch-and-bound method can be stopped when the first 
feasible solution for P is found. In this way a good 
solution, but not necessarily the best one. is obtained. In 
most applications this is no problem because the coefficients 
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in the model are estimates and the differences between the 
exact solution and the one found can be made arbitrarily- 
small . 

The heuristic also uses the reduced costs to fix 
nonbasic variables to 1 or 0: 

1) Fix xj to 0 if xj = 0 in the continuous solution and 

Z LP - K 2 Z LP < d j 

2) Fix Xj to 1 if Xj s 1 in the continuous solution and 

Z LP - K 2 Z LP < -dj 
where K2 < 1. The above rules are applied after the 

continuous solution of the relaxation of P is found. The 
value cf cannot be chosen as high as the value of K2. 
because when choosing it must be certain that the value of 
the objective function for the solution of P is larger than 
^1 Z LP- ^ the value of or K2 is too large, the decision 
tree is small. Then it does not take much time before it is 
clear that no solution to P can be found for the chosen 
values of and * n suc h a case the values of and/or 

K2 can be adjusted and the procedure is started anew. The 
following example illustrates the efficiency of the above 
heuristic . 
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Example 

In achievement test construction a high value for the 
reliability coefficient of the test is wanted and this goal 
can be achieved by selecting items from a test item bank with 
large contribution to the reliability (van der Linden & 

9 
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Adema, 1987). However, other test construction goals are also 
possible (Theunissen, 1985; van der Linden & Boekkooi- 
Timminga, 1988). In practice, test item banks usually consist 
of hundreds of items and practical constraints have to be 
imposed on selection of items. 

The heuristic has been applied to zero-one programming 
models with the goal as mentioned above in the objective 
function. The (0 , 1 )-variables x-j were defined as follow: 



xi = i 



0 item j not in the test 

1 item j in the test. 



The total time of test administration and the mean of the 
difficulties of the items in the test were restricted. Except 
for these constraints all the coefficients in the constraints 
were -1, 0 or 1 . One constraint was introduced to fix the 
number of items in the test Also the numbers of items to be 
selected from different subdomains of the item bank were 
restricted. Sometimes we want two items to be simultaneously 
included in or excluded from the test. Constraints to satisfy 
this wish were included in the model. In all, 14 constraints 
were imposed. 

A number of test construction problems were implemented 
and solved by the proposed heuristic for item banks with 400 
items. This was done on a DEC2060 computer. The modifications 
were introduced in the program LANDO. It was assumed that the 
items in the item bank satisfy an item response modex. An 
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item response model specifies a relationship between the 
observable examinee test performance and the unobservable 
trait or ability assumed to underlie performance on the test. 
The relationship between the "observable" and the 
"unobservable'* quantities is described by a mathematical 
function. For this reason, item response models are 
mathematical models based on assumptions about the test data 
(Hambleton & Swaminathan, 1985). In the implementations, the 
Rasch model and the 3-parameter logistic model were used. The 
probability that an item i is answered correctly by a person 
with ability 0 under the Rasch model is 

-1 

P(+| i,0) = [1 + exp (b i - 9)] 

where is the difficulty of item i and repesents the point 
on the ability scale at which an examinee has a 50 percent 
probability of answering item i correctly. Under the 3- 
parameter model this probability is 

P(+|i.e) = c i + id - CjL ) I (1 + exp (-a^G -b^)] 

where the parameter c^ represents the probability of 
examinees with low ability correctly answering an item and 
the parameter aj_, called item discrimination, is proportional 
to the slope of P(+|i,9) at the point 9 = b±. 

The results are shown in Table 1. To see how important 
the choice of is, two values, .99 and .995, were chosen 
for . Parameter K2 was set equal to .998 if 20 items were 

ERLC 1 1 
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selected and equal to .999 if 40 items were selected. It is 
possible to choose K2 higher when 40 items are selected, 
because the difference between Z + and Z^p in percents is 
smaller for 40 items selected than for 20 item selected. 



Insert Table 1 here 



For the 3-parameter model the dispersion in the 
coefficients of the objective function was greater. 
Therefore, more variables were fixed after applying the rules 
with the reduced costs. As a conseqence, the CPU-times for 
the 3-paraineter model were better. 

Because of the modification in the initialization, more 
branches are fathomed without finding a feasible solution o£ 
P for Kj = .995 then for Kj = .99. Therefore it is possible 
that more CPU-time is needed to solve P for = .995. 

Conclusions 

A heuristic for solving large-scale zero-one programming 
problems is proposed. This heuristic is useful in particular 
when the optimal objective function value of P is close to 
Z^p. The heuristic was used to solve test construction 
problems. As shown in Table 1, the CPU-times for solving the 
zero-one programming problem were close to the CPU-times 
needed to solve the continuous problem. 

ERLC * 2 
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Table 1 

CPU-times and objective function values for different values 
of and K2 









Zero— one 


Continuous 


ni 


1 




Obiective CPU 


Objective 


CPU 








Function (sec) 


Function 


(sec ) 4 








Value 


Value 










Rasch model 






20 


. 990 


. 998 


8 80 22 17 


8.8533 


18.04 


20 


.995 


.998 


8 82 24 65 


8.8533 


18.04 


40 


.990 


.999 


17.27 34.82 


17.3911 


24.63 


40 


.995 


.999 


17.32 35.96 


17.3911 


24.63 








3-parameter model 






20 


.990 


.998 


9.43 18.35 


9.4468 


16.68 


20 


.995 


.998 


9.43 18.93 


9.4468 


16.68 


40 


.990 


.999 


18.35 25.46 


18.3500 


25.42 


40 


.995 


.999 


18.35 25.47 


18.3500 


25.42 



Note , ni = number of items in the test 
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